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Despite deleterious effects on individuals, the t haplotype is a selfish genetic element present in many house mouse populations. 
By distorting the transmission ratio, +/t males transmit the t haplotype to up to 90% of their offspring. However, t/t individuals 
perish in utero. Theoretical models based on these properties predict a much higher t frequency than observed, leading to the t 
paradox. Here, we use empirical field data and theoretical approaches to investigate whether polyandry is a female counterstrategy 
against the negative fitness consequences of such distorters. We found a significant decrease of the f frequency over a period of 
5.5 years that cannot be explained by the effect of transmission ratio distortion and recessive lethals, despite significantly higher 
life expectancy of +/t females compared to +/+ females. We quantified life-history data and homozygous and heterozygous 
fitness effects. Population subdivision and inbreeding were excluded as evolutionary forces influencing the t system. The possible 
influence of polyandry on the t system was then investigated by applying a stochastic model to this situation. Simulations show 
that polyandry can explain the observed t dynamics, making it a biologically plausible explanation for low t frequencies in natural 
populations in general. 
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In its classical conception, Darwinian evolution by natural se- 
lection predicts that genes have to contribute to organismal 
fitness to be successful. An increasing number of cases are emerg- 
ing where this paradigm is violated. Selfish genetic elements 
define such heritable entities. They spread through populations 
despite being associated with negative fitness consequences for 
the organism (Burt and Trivers 2006). These stretches of DNA 
distort Mendelian segregation in their favor (transmission ratio 
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distortion (TRD) or meiotic drive) and thereby gain an advantage 
over their wild-type variants. Selfish genetic elements fascinate 
evolutionary biologists because they demonstrate persuasively, 
that systematic advantages at a gene or gamete level can lead to a 
distinct disadvantage at a higher level of organization such as the 
individual or population (Okasha 2006). 

The t haplotype. Since its discovery in 1927 
(Dobrovolskaia-Zavadskaia and Kobozieff 1927), the t haplotype 
in house mice (Mus domesticus) has been intensively studied and 
is now the best known example of transmission ratio distorters. 
Efforts to understand the underlying genetics of this distortion 
have been quite successful. Heterozygote males produce equal 
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proportions of + and t sperm (Silver and Olds-Clarke 1984), but 
an elaborate poison antidote system exclusively impairs flagel- 
lar function of the + sperm (Schimenti 2000; Lyon 2003; Bauer 
et al. 2005, 2007). This leads to a transmission of t gametes 
that deviates considerably from 50%, and up to as much as 99% 
(Silver 1993; Lyon 2003). Females, on the other hand, transmit 
the t gametes in the usual Mendelian ratios. The t haplotype com- 
prises a complex of linked genes as large as 20 cM (30^-0 Mbp), 
occupying the proximal third of chromosome 17 (Silver 1993). 
Four major, nonoverlapping inversions block recombination and 
assure that the complex is transmitted as one intact entity (Artzt 
et al. 1982). It has been found in all four house mouse subspecies 
across the world and is thought to have existed in house mouse 
populations for 1.5 to 2 million years (Silver 1993). As is the 
case for most known distorter systems, they have negative fitness 
effects on an individual level. Cases with positive or neutral ef- 
fects are probably not detected because they will drive to fixation 
in very short time. Numerous variants of t haplotypes have been 
found. Most t haplotypes carry lethal recessive mutations such 
that homozygotes for the same variant perish in utero. On the 
other hand, homozygosity for different, complementing variants 
(e.g., t x /t y ) invariably result in male sterility (Klein et al. 1984). 

The t frequency paradox. Despite the solid understand- 
ing of the structure and mechanisms, the implications for popu- 
lation dynamics of this genetic polymorphism are still a puzzle 
(see Ardlie (1998) for a review). In one of the first theoretical 
approaches to this system more than 50 years ago, Bruck (1957) 
provided a mathematical model that took into account the an- 
tagonism between segregation distortion in males supporting the 
t haplotype frequency, and viability selection at the embryonic 
stage acting against it. The model assumes random mating and 
an infinite, unstructured population. Under these assumptions, the 
frequency of the t haplotype reaches a steady-state equilibrium 
that is surprisingly high, despite the strong negative fitness conse- 
quences for t/t individuals. With a TRD of 0.9, an equilibrium t 
frequency of 0.33 is predicted, meaning that two-thirds of a pop- 
ulation are expected to be +/t heterozygotes. This expectation is 
however in marked contrast to frequencies usually found in wild 
populations. Various empirical studies measuring the frequency of 
t haplotypes in natural populations of different subspecies around 
the globe show the same picture: frequencies are on a persistent 
but low level, with t frequencies ranging between 0.05 and 0.15 
(Ardlie and Silver 1998; Lenington et al. 1988; Huang et al. 2001; 
Dod et al. 2003). This discrepancy — often called the t frequency 
paradox — has evoked a large number of theoretical models (both 
analytical models and stochastic simulations) focusing on vari- 
ous forces that might account for this low frequency in a natural 
context. 

Population subdivision. Population subdivision and ge- 
netic drift were considered as possible factors to reduce t 



haplotype frequency, as they reduce heterozygosity and thereby 
the frequency of t carrying individuals (Lewontin and Dunn 1 960; 
Levin et al. 1969; Petras and Topping 1983; Nunney 1993). It has 
nevertheless been a matter of debate whether mouse populations 
are as substructured as assumed by these models. Some studies 
point out that overall t frequencies within a population are rather 
stable (Lenington et al. 1988), whereas empirical data indicate 
that population size has an influence on t allele frequency, with 
frequencies decreasing in larger populations (Ardlie and Silver 
1998). 

Heterozygote fitness effects. Selection against + /t het- 
erozygotes will substantially reduce t frequencies (Young 1967; 
Lewontin 1968; Johnston and Brown 1969; Hartl 1970). Here 
again, empirical data do not show a clear picture. Besides the 
expected reduction in litter size after a double +/t mating, sev- 
eral studies also showed litter size reductions if only one sex was 
carrying a t (Johnston and Brown 1969; Lenington et al. 1994; 
Carroll et al. 2004). Viability, on the other hand, has been found 
to be either higher (Dunn et al. 1958) or lower (Carroll et al. 2004) 
in +/t heterozygotes independent of sex. 

Female choice. In a series of studies, Lenington and col- 
laborators found evidence for sexual selection against t carry- 
ing males via female preferences (Lenington and Egid 1985; 
Lenington et al. 1994). Using a Y-maze apparatus, they showed 
repeatedly that +/t females spend significantly more time near 
+/+ male-derived odor cues when given the choice between +/+ 
and +/t. This preference is adaptive in light of the reduction in 
litter size that would result from a mating with a +/t male. How- 
ever, the question whether these odor preferences reflect actual 
mating preference is unresolved. Furthermore, male dominance 
status seems to be an even more important predictor for female ol- 
factory preference (Coopersmith and Lenington 1992), although 
evidence of the effect of t on male dominance is contradictory 
(Lenington et al. 1996; Carroll et al. 2004). 

Genetic modifiers. Theoretical work predicts the evo- 
lution of genetic modifiers suppressing segregation distortion 
(Charlesworth and Hartl 1978). Such modifiers have been found 
in various other distorter systems (Hiraizumi and Thomas 1984; 
Atlan et al. 2003). Surprisingly, they appear to be rare for / hap- 
lotypes of natural house mouse populations (Ardlie and Silver 
1996), although there are reports on suppressors from laboratory 
mice (Bennett et al. 1983; Gummere et al. 1986). 

Polyandry. Even though a reduction of the segregation ra- 
tio due to modifiers seems to be nearly absent in wild populations, 
there still remains the possibility of a reduction of TRD by other 
means, such as polyandry. Haig and Bergstrom (1995) proposed 
the idea that females could avoid individual fitness reduction 
due to distortion by systematic multiple mating. The hypothe- 
sis of polyandry as a general counterstrategy against distorters is 
based on the observation that gametes carrying the t haplotype are 
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"by definition" strong intraejaculate competitors, but do compar- 
atively poorly in competition with other ejaculates (Zeh and Zeh 
1996; Zeh and Zeh 1997; Price and Wedell 2008). In comparison 
to previously described precopulatory mating preferences, this 
mechanism would not even require the female ability to distin- 
guish between the two different genotypes. Hence, polyandry is a 
simpler and potentially more robust strategy to counter the spread 
of t haplotypes. A series of empirical studies on several species, 
predominantly on the genus Drosophila, support this idea. Both 
fertility reductions and reduced competitive ability of sperm from 
males carrying selfish genetic elements have been found in var- 
ious systems (see Price and Wedell (2008) for an overview). In 
Drosophila pseudoobscura the X-linked driver (SR) was found to 
reduce sperm competitive ability (Price et al. 2008a) and females 
evolved increased remating rates in the presence of the distorter 
(Price et al. 2008b). Fertility reduction and negative effects on 
sperm competitive abilities were also found in the SR system and 
Wolbachia of Drosophila simulans (Snook et al. 2000; Atlan et al. 
2004; Champion de Crespigny and Wedell 2006). 

This hypothesis appears to be quite promising for the house 
mouse case, because recent studies indicate that female house 
mice are indeed actively polyandrous both in the wild (Dean 
et al. 2006) and in an experimental context (Rolland et al. 2003). 
Polyandry provides opportunity for postcopulatory selection pro- 
cesses such as sperm competition and cryptic female choice. It has 
already been shown that polyandry increases offspring postbirth 
survival (Firman and Simmons 2008c) and facilitates inbreeding 
avoidance (Firman and Simmons 2008b). As mentioned above, 
evidence for negative effects of the t on male fertility are con- 
vincing: apart from t/t sterility (in the case of complementing 
t variants), reduced fertility (of about 20%) has invariably been 
reported in +/t heterozygous males (Johnston and Brown 1969; 
Lenington et al. 1994; Carroll et al. 2004). However, empirical 
data on sperm competitive abilities in relation to the t haplotype 
are still scarce. One study looking at paternities in a wild popula- 
tion yielded a mean fraction of 0.17 +/t among litters involving 
both +/t and +/+ fathers (based on three litters, Ardlie and 
Silver (1996)), another study using controlled sperm mixing ex- 
periments obtained a +/t proportion of 0.22 (based on eight litters, 
Olds-Clarke and Peitz (1986)). In the study of Carroll et al. (2004) 
on seminatural enclosure populations, the t was transmitted to 
36% of the offspring. 

Aims. The present study focuses on both empirical and the- 
oretical approaches to understanding the dynamics of the t haplo- 
type on a specific wild house mouse population. First, we provide 
a theoretical model looking at possible impacts of polyandry on 
t haplotype frequency. To our knowledge, no theoretical model 
has investigated this for t haplotypes. However, the evidence in 
the previous section make polyandry a promising evolutionary 
force to explain low t frequencies in natural populations. Second, 



extensive data collection in a free-living population of house mice 
over a time period of more than five years allows us to estimate 
many parameters likely to affect the t frequency. In addition to 
well-described parameters such as distortion level and homozy- 
gous fitness effects, we were able to obtain reliable estimates on 
heterozygous fitness effects, the degree of inbreeding, as well as 
t frequency dynamics. Reliable estimates of these have not been 
available for natural house mouse populations so far (Burt and 
Trivers 2006). In addition, we were able to describe a series of 
general life-history parameters (generation time, life expectancy, 
net reproductive rate) as yet largely unknown for wild house mice. 
In a third step, the theoretical model and the empirical data were 
compared. We ran computer simulations using the parameters es- 
timated from our population to predict the t frequency dynamics. 
With the present approach, we are able to test specific models with 
parameters directly estimated from the population of interest. 

Materials and Methods 

THE MODEL 

A classical Fisher- Wright population with infinite population size 
and no mutation is assumed. For the purposes of the model, the 
whole t haplotype is simplified to a single locus with two alleles — 
the wild-type allele + and the distorter allele t. Using i = 1 for 
females and i = 2 for males, the variables and for 
each sex denote the frequency of alleles + and f, respectively. 
Similarly, for each sex i, the variables P++(i), P+t(i), and P ff(i ) 
describe the frequency of the genotypes +/+, +/f, and t/t, re- 
spectively. For example P++(\) would represent the genotype fre- 
quency of female homozygote wild-types, t/t homozygotes are 
inviable and it simply holds that p t( j) = | P+t(i) and P ++(,■) = 1 — 
P+t(i)- This also explains why p«j) can never exceed 0.5, the case 
in which all members of the population are heterozygotes. The 
sex-independent t allele frequency p t can be expressed as p t = 
fiPt(i) + fiPt(2), where f\ and ft are relative frequencies of 
females and males in the population (f i + f 2 = 1). 

The Life Cycle. To determine the allele frequencies p',^ in 
the next, nonoverlapping generation, the following life cycle was 
used. 

(1) First, adult individuals of the present generation mate 
with each other. Mating probabilities are based on the 
genotype frequencies, which means that all matings are 
random. A certain fraction of the females 0 < \|x ^ 1 is 
assumed to mate twice, whereas the rest of the female 
population 1 — v|/ mate with one randomly encoun- 
tered male. i|/ is independent of female genotype. It is 
also assumed that matings are not male limited, that 
is, the number of males does not influence the mating 
frequencies. The four monotypic matings (two female 
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genotypes can each mate with two male genotypes) and 
six additional polyandrous matings sum to 10 possible 
mating combinations (see Table 2). Mating frequencies, 
expressed as probabilities, can be calculated by multi- 
plying the parental genotype frequencies. For example, 
the frequency of single matings between a female of 
genotype a and a male of genotype b can be expressed 
as 

p(d) n p(b) = (l - Wd)^), (i) 

where a, b G {+ + , + t }. 

The frequency for a multiple mating involving the or- 
dered male genotypes b\ and Z?2 are given by 

P{a) n P(bl) n P(b 2 ) = tyPa { l)P bl (2)Pb 2 <2)- (2) 

Note that equations ( 1 ) and (2) determine the probability 
for all possible matings, but they do not determine the 
outcome of the matings. This is done in the next section. 
(2) After mating, the proportion of the different zygotes of 
the different crosses and double-crosses is calculated. 
We distinguish between within and between ejaculate 
effects. 

Within ejaculate effects. Because meiotic drive changes 
the proportion of functional gametes in heterozygous 
males, the expected zygotes of the affected crosses devi- 
ate from Mendelian predictions. Segregation ratio (pro- 
portion of t gametes in the functional male gamete pool) 
in males is characterized by the variable 0 < x < 1 . 
Between ejaculates effects. If a female mates with more 
than one male, both "quantity" and "quality" of the 
sperm are assumed to determine the fertilization suc- 
cess of the involved males (gamete fitness). The mating 
order of the males shall not have an influence on fertil- 
ization success. Quantity reduction, represented by the 
coefficient v, results from the reduction of functional 
gametes through ratio distortion x. Under the assump- 
tion that every male produces the same amount of sperm 
(functional and dysfunctional) per ejaculate, the fraction 
of functional sperm in a given ejaculate is dependent on 
TRD because the latter operates on rendering + sperm 
dysfunctional. In +/+ males that do not exhibit TRD, 
no sperm is lost and thus v ++ = 1. In +/t males, v is 
dependent on the level of x and can be expressed by 
v +t (x) = (1 + \2x — 1|) . The higher the deviation 
from x = 0.5, the bigger the individual males' loss in 
functional sperm. For x > 0.5, it can also be described as 
v +t (x) = l/2x. Because this expression summarizes the 
relevant cases for the t haplotype and is mathematically 
easier to handle, it will be used here. 



Differences in quality between the two males ejaculates 
are modeled with the parameter c, which describes the 
relative disadvantage of the remaining heterozygote de- 
rived sperm. Hence, relative competitiveness of +/+ 
male sperm is w ++ = 1, relative competitiveness of 
+/t male sperm w +t = 1 — c. 

Given two males of genotype b\, b 2 G {+ + , + t }, 
the probability F bl of fertilization of a given egg by the 
male of genotype b\ is dependent on both gamete fitness 
components w and v and can be expressed by 



w bl v bl + w bl v bl 

If both males have the same genotype (e.g., b\=b 2 = 
'+ f ')> men Fh = 0.5. Note that F bl + F bl = 1. Know- 
ing this, the proportions of zygote genotypes produced 
by each mating cross P~ can be determined (see Table 1 
for one illustrative example and Table 2 for an overview 
over all possible mating crosses). The zygote frequen- 
cies of the total population PL^ of the genotypes a G 
{++, + ?, tt} are then given by the sum of the out- 
comes of each individual mating cross P~~ weighed by 
the frequency of the respective mating cross derived in 
equations (1) and (2). 
(3) In the next step, viability selection occurs. The different 
genotypes are given different probabilities to establish 
themselves as adults in the population. Estimates from 
our wild population suggest that these probabilities can 
be sex dependent and eventually lead to sexually antag- 
onistic effects. Relative viability differences between 
+/+ and +/t individuals will be characterized by s ; , 
where i again defines sex. Overdominance is indicated 

Table 1. One example showing how to calculate the zygote fre- 
quencies Pi of genotypes a e {+ +, + t, tt} for the multiple cross 
involving a female of genotype a = + f and males of genotype 
b\ = ++ and bi = + f. F^ and F^ are derived from equation 
(3). The zygote frequencies for all other crosses can be calculated 
likewise. 



Male b\ 


Male b 2 




Gametes + 


+ 


t Total 


Female a + ^F bl 


\F bl (\-x) 


l -F, x ± 

2 *2 2 


t \F bl 


\F bl (\-x) 


if, f I 
2 1 h 1 2 


Total F bl 


Fbi 


1 


Resulting zygote frequencies ( 


using eq. 3) 




P+/+ = 2 F b, + 2~F hl (\ ~X) = 


l— c+x+xc 
2(l-c+2t) 




Pl/t = iFbi + \F bl {\ — x)+l 


F bl x = \ 




pz l 77 T x—xc 
r t/r — 2 r '! L ~ 2(1-c+2t) 
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Table 2. The different possible single (female of genotype a mating with a male of genotype b) and multiple crosses (female of genotype 
a mating with males of genotypes bi, bz) in the polyandry model, including their resulting zygote proportions (P* for a e {+ +, + t, 
tt}). Expected paternity shares at birth (correcting for t/t lethality) of the first male by, and the second male bi, are given in the cases of 
multiple crosses. 



Zygote frequencies per mating cross Paternity shares at birth 







p- 7 . . 
r ++ 


1 +t 


p : 




"2 




+/t * +/t 


l— t 

2 


1 

2 


1/2 






Single matings 


+/+ * +/t 




1 

2 


0 






male b x female a 


+/t * +/+ 


1 — 1 


1 


0 








+/+ * +/+ 




0 


0 








(+/t, +/t) * +/t 


1— t 

2 


1 

2 


1/2 


1 

2 


l 

2 




(+/+, +/+) * +/t 




1 

2 


0 


1 

2 


1 
2 


Multiple matings 


(+/+, +/t) * +/t 


1— C+T + TC 


1 


x— tc 


4x 


2-T-2c+Tf 


2(l-c+2t) 


2 


2(l-c+2t) 


2+3t— 2c+xc 


2+3t-2c+tc 


(male b\, male b-i) x female a 


(+/+, +/+) * +/+ 




0 


0 


1 

2 


1 
2 




(+/t, +/t) * +/+ 


1 — 1 


1 


0 


1 

2 


1 

2 




(+/+, +/t) * +/+ 


1— c+x+xc 


x—xc 


0 


2t 


1-c 






1 _^-l-9 T 


1 







by < 0 and underdominance by s t > 0. Thus, rela- 
tive viabilities for +/+ genotypes are w ++ ^ = 1 and 
w + ,d) = 1 — Si for +/t genotypes, respectively. Because 
the t haplotypes of our study population carry an em- 
bryonic lethal (A. K. Lindholm, unpubl. data), viability 
of homozygote t carriers is set to be Wu(f) = 0. Adult 
genotype frequencies for any genotype a &{+ + , + t} 
after viability selection are given by 



We use a scenario in which \|/ = 0, c = 0, and s,- = 
0 as a null model. This is equivalent to the model by 
Bruck (1957), and is a case of dominance of the wild- 
type allele, because both +/+ and +/t have the same 
fitness, whereas t/t are lethal. 

THE STUDY POPULATION 

All the data used in the present study originate from a free-living 
population inhabiting a barn near Zurich. The population was es- 
tablished in 2002 from 12 founder individuals (originating from 
two different capture sites close-by). On an area of 72 {m 2 , an- 
imals were given breeding opportunities (40 nest boxes) and ad 
libitum food and water (around eight drinking and feeding sites). 
Branches and shelves allowed mice to establish several territories 
in the building and freely enter and leave the barn through nu- 
merous openings. The population has been intensively monitored 
from inception for births and deaths, and tissue samples for ge- 
netical analyses have been taken from each pup and adult, as well 
as from individuals found dead. For the period between 2003 and 



summer 2008 used for the present study, 2177 pups were sampled 
at an age of 13-15 days. We found a rate of about 5% unsampled 
individuals among all captured adults during this study period. 
We thus sampled at least 95% of all pups born in the population. 

GENOTYPING 

To identify the t haplotype, the Hba-ps4 locus — a marker contain- 
ing a 16-bp t haplotype-specific insertion (Hammer et al. 1989) — 
was amplified and scored. Sexing of pups was performed by 
amplification of three different Y-specific microsatellite mark- 
ers (Y8, Y12, and Y21; Hardouin et al. (2010) and M. Teschke, 
pers. comm.). Samples in which no Y marker amplified were 
scored as females (A. K. Lindholm, unpubl. data). Identification 
of dead and adult individuals was also carried out by matching 
multi-locus genotypes of the dead or adult individuals with those 
of the sampled pups at 21 unlinked microsatellite loci, allowing 
one mismatch (A. K. Lindholm, unpubl. data). As pups were sam- 
pled at an age of 13 days, birth date could be calculated. Using the 
date on which the dead body was found as date of death, longevity 
could be determined. 

Results 

PREDICTIONS BASED ON DETERMINISTIC MODEL 

To investigate the potential impact of polyandry on the t frequency 
p,, the life cycle described by equations (1-4) was repeated for 
1000 generations. Independent of the initial parameter settings, 
p t reached a steady-state equilibrium p t after only a few genera- 
tions. Figure 1 shows p t dependent on multiple mating frequency 
\|/ and heterozygote relative sperm competitiveness c. Segrega- 
tion distortion was set to the level observed in our population at 
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0 0 0 2 0 4 0 6 0 8 1 0 

multiple mating frequency 



Figure 1. Equilibrium t frequency pt dependent on multiple mat- 
ing frequency (i|r) for different levels of t sperm competitiveness 
-1 ^ c < 1 and x = 0.90. The solid line represents the case 
where heterozygote males only suffer from quantitative sperm 
disadvantages resulting directly from segregation ratio distortion 
(c=0). 

x = 0.90 (see Supporting infor). If +/t males have only a quanti- 
tative sperm competition disadvantage (lower viable sperm count) 
compared to +/+ males (c = 0, solid line in Fig. 1), polyandry 
does not have strong effects on t frequency. When all females mate 
multiply (ijr = 1), the proportion of +/t individuals in a population 
is reduced by about 20% compared to the case where only single 
matings occur (\[r = 0). Frequency of matings involving a +/t and 
a +/+ male — the cases where sperm competition disadvantage 
can play a role at all — is obviously too rare to produce substantial 
effects on p t . This changes however, if disadvantages in sperm 
quality are included in the quantitative sperm competition disad- 
vantage of t carrying males (c > 0) . Higher disadvantages in +/t 
sperm quality lead to lower equilibrium t frequencies. Depend- 
ing on the level of c, effects on t frequency can be quite strong, 
especially when accompanied by high multiple mating rates. 

PREDICTIONS BASED ON STOCHASTIC MODEL AND 
COMPARISON TO WILD HOUSE MOUSE POPULATION 

Empirical observations. The model was applied to an inten- 
sively studied, free-living house mouse population near Zurich. 
A lethal version of the t haplotype was present in this population 
since its establishment. Figure 2A shows the frequency of the 
t over a 5.5-year time period (2003 until summer 2008) among 
2177 pups. Although fluctuating substantially, t frequency de- 
creased significantly over these five years (GLM on proportions 
using a binomial error structure and a logit link function using 



time as a continuous factor (in years); residual deviance: Xn = 
79.53, P < 0.001). 

The t paradox. This decrease is surprising given the sub- 
stantial advantage for the t haplotype through segregation ratio 
distortion. For the levels of segregation ratio distortion found in 
the population (x = 0.90, see Supporting information for more de- 
tails), Bruck (1957) predicts an equilibrium t frequency p, as high 
as 0.32 (see Fig. 2A). In addition, we found that t heterozygote 
females have a significantly longer life expectancy than homozy- 
gote wild-type females (Cox proportional hazard model: n = 174, 
exp(f5)=2.46, P < 0.01, see Figure 3 and Supporting information 
for more details). Male survival on the other hand was not depen- 
dent on the t haplotype (Cox proportional hazard model: n = 185, 
exp(f5)=1.26, P = 0.31). The higher survival probabilities of fe- 
male heterozygotes hence leads to overdominance (as compared 
to the Null model, which is a case of dominance). We used data 
from 21 unlinked microsatellite markers to examine inbreeding. 
Inbreeding levels did not change throughout the whole obser- 
vation period. A linear regression analysis showed that the F IS 
(deviations in heterozygosity from Hardy- Weinberg predictions) 
did not significantly change over time (Regression coefficients: 
a = 0.0259, p = -0.0065, n = 147, P = 0.07, r = 0.023, see 
Fig. SI). 

Polyandry. Based on paternity analysis, we find a propor- 
tion of litters sired by more than one male (multiple paternity 
rate) of around 0.3 (Camani 2005). Competitive ability c of +/t 
derived sperm was estimated from laboratory experiments with 
wild mice (A. Manser, unpubl. data). Using an elaborate choice 
device that allowed females free physical access to both a +/+ 
and a +// male without male interference, we found a propor- 
tion of 0.19 (95% CI using a binomial error distribution: [0.08, 
0.35]) sired by the +/t male among four multiply sired litters. The 
experimental setup did not control for mating order. This would 
correspond to c = 0.58 (following Table 2). We are, however, 
facing a fundamental problem. Multiple paternity rate and sperm 
competitiveness of c = 0.58 mark only the "minimal" amount of 
sperm competition possible in our system. Cases of multiple mat- 
ings with only one successful competitor do not lead to multiple 
paternity and are therefore missed by our analysis. The amount of 
such cases highly depends on sperm competitiveness c, the litter 
size X as well as t frequency p, in the population. But we can- 
not estimate both c and \|/ only knowing multiple paternity rate. 
We are left with an underdefined system that forces us to make 
assumptions. 

Scenario I: This most conservative scenario assumes that all 
multiple matings resulted in multiple paternities. 
This leaves us with a multiple mating frequency 
of \[r = 0.3 and a t sperm competitiveness of 
c = 0.58. 
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Null Model: dominance, no polyandry (v=0. c=0) 
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Figure 2. Observations from the wild population and simulation predictions from October 2002 until June 2008. The first dot within 
vertical dashed lines represents the initial t frequency among the first 50 adults ± SE. The following dots represent the f frequency among 
the pups born in the given 3 month time interval ± SE (n (ot — 2177). Solid lines show model predictions with 95%-confidence intervals 
(shadows) for (A) the null model (dominance and no polyandry) and (B)-(D) scenarios l-lll (in Table 3, overdominance and different levels 
of polyandry). [Correction added June 1, 2011 after Online publication: Figure 2 legend updated to reflect black-and-white publication.] 
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Figure 3. Empirically estimated survival functions for (A) females and (B) males of different genetic backgrounds based on a Cox model. 
Each estimate is accompanied by a pointwise 95% confidence envelope. 



Scenario II: In this intermediate scenario, multiple mating fre- 
quency is assumed to be slightly higher. v|/ = 0.6 is a 
value that is within the range of what has previously 
been described for wild house mouse populations 
(Dean et al. 2006). By making use of our model 
(with known litter size X and mean t frequency p t ) 
we find that c = 0.88 is needed to explain a gap 
between multiple paternity and multiple mating of 
this order. 

Scenario III: The third scenario assumes an almost maximal mul- 
tiple mating rate. It is based on the only study that 
actually looked at the mating behavior of females in 
an experimental context (Rolland et al. 2003). Given 
the choice between two different males, they found 
that 20 of 21 females (95%) mated with both of the 
males. Therefore, this scenario assumes \|/ = 0.95 
and c = 1 . 

Stochastic model. To determine if polyandry can explain 
this dramatic decrease in frequency, we ran simulations based on 
parameters estimated from the population (see Table 3 and Sup- 
porting information for a more detailed description). The deter- 
ministic, panmictic model described before was slightly modified. 
The deterministic nature of the model assuming infinite popula- 
tion size was now replaced by a stochastic simulation sampling 
randomly through the different life cycle stages based on effective 
population size and mean litter size. Life-history data from the 
populations allowed us to calculate an average time to reproduc- 
tion (generation time) of 9 months (see Supporting information). 
This estimate was used to fit the simulation assuming nonover- 
lapping generations to the study population, where generations 
overlap. To get reliable estimates on model predictions and its 



Table 3. Parameter estimates used in this study. Roman numerals 
refer to the three scenarios tested. Estimation of the lower part is 
shown in the Suppo. 





Definitions 


Estimates 




Frequency of females mating twice 


I: 0.30 


II: 0.60 


III: 0.95 


c 


Relative disadvantage of +/t derived 
sperm 


I: 0.58 


II: 0.88 


III: 1.00 


X 


Average litter size at birth 


5.47 


Y 


Generation time 


9 months 


N e 


Effective population size 


50 


X 


Segregation ratio distortion 


0.90 


Sl 


Difference in survival between + /+ 
females and +/t females 


-0.30 


S2 


Difference in survival between +/+ 
males and +/t males 


0.00 



confidence intervals, the 5.5-year observational period was sim- 
ulated 10,000 times for each scenario. Initial frequency was set 
to the t frequency among the first 50 adults (equal to N e , shown 
in red in Fig. 2, separated by a dotted line) which was p, = 
0.24. 

Figure 2 A shows the stochastic model predictions for the null 
model (following Bruck (1957)), which assumes no polyandry, 
equal fitness between +/+ and +/t individuals, and homozygote 
t/t lethality (dominance). Figure 2B-D shows the three differ- 
ent polyandry scenarios that include female overdominance and 
sperm competition. Scenarios 2 and 3 with high polyandry fit the 
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data considerably better than the null model. Notice that — in con- 
trast to scenario 2 — the deterministic equilibrium prediction for 
scenario 3 is the extinction of the t. However, more generations 
than used in the present simulation are on average needed to reach 
extinction. 

EFFECTS ON MEAN POPULATION FITNESS 

Examining mean fitness of the overall population here reveals 
interesting effects of the t haplotype. As is the case for all self- 
ish genetic elements (Burt and Trivers 2006), the t haplotype 
elicits conflict between different levels of selection (allelic prop- 
agation versus population fitness). This conflict can be revealed 
when looking at the mean population fitness and its related gene 
frequencies. In the case of overdominance, mean population fit- 
ness is generally optimized by maximizing the number of het- 
erozygotes in a population (Wright 1929). However, in the pres- 
ence of distorters, the expectation that Mendelian segregation and 
random mating lead to gene frequency changes that increase mean 
population fitness (Wright 1929) no longer holds. By increasing 
their own frequency, t haplotypes drag the population frequency 
away from the optimum. This is shown in Figure 4 for our specific 
case: in accordance with Wright, without segregation distortion 
(solid line) an equilibrium frequency that optimizes population 
fitness is reached (point [1]). As soon as we add meiotic drive 
to the system (dashed line), mean population fitness at equilib- 
rium is far from the optimum (point [2]), manifesting the conflict 
between gene and population level fitness. Dependent on the sce- 
nario, these population fitness losses are recovered by the process 
of sperm competition. Whereas scenario 2 has an equilibrium t 
frequency that nearly corresponds to the optimum of population 
fitness, the deterministic model predicts extinction of the / and 
suboptimal population fitness in scenario 3. The observed t fre- 
quency at the end of the observation period (shown as vertical 
dash-dotted line) falls between maximum population fitness and 
extinction of the t haplotype. 

Discussion 

t dynamics. One of the main findings of this study is the signifi- 
cant decrease in t frequency in a wild population over a time period 
of five years. It is remarkable how strong the dynamics of the t 
haplotype in our study population closely resemble the patterns 
observed in other studies that examined t frequencies in a natural 
or seminatural context. Both Ardlie and Silver (1998) and Carroll 
et al. (2004) find low frequencies in their study populations of 
house mice, the latter report a similar decrease in frequency over 
time in their enclosure populations. Even in a different subspecies 
(Mus castaneus) in Taiwan, t frequencies are found at persistent 
but low frequencies (Huang et al. 2001). These striking parallels 
among studies measuring t frequencies in different geographic 



and phylogenetic contexts suggest similar general mechanisms. 
It also confirms our assumption that if we are able to unravel 
the factors determining the t dynamics in our specific population, 
these findings have implications for house mouse populations in 
general. 

Parallels to other studies are not only restricted to t frequency 
dynamics, as most other parameter estimates confirm what has 
previously been found. TRD of 0.90 is a level typically observed 
in other natural and semi-natural populations (Ardlie and Silver 
1996; Huang et al. 2001; Carroll et al. 2004), and provides no 
indication of a reduction of distortion level by modifiers (Ardlie 
1998). The estimates concerning polyandry in our study popula- 
tion are also consistent with previous investigations. A marked 
multiple paternity rate of about 30% (Camani 2005) provides 
further support that house mice mating systems are highly polyg- 
ynandrous. This estimate is further in high accordance with recent 
work, that found multiple paternity rates of 23% and 26% (Dean 
et al. 2006; Firman and Simmons 2008a). 

Competitiveness of the +/t derived sperm c is clearly the 
model parameter with the weakest support, both in our study and 
in the literature. Two other studies measured +/t male paternity 
shares of 0.17 and 0.22 (Olds-Clarke and Peitz 1986; Ardlie and 
Silver 1996), proportions that are comparable to what has been 
found here (0.19). They contradict however slightly with Carroll 
et al. (2004) who found a paternity share for +/t males of 0.36, 
a value that would suggest no "qualitative" differences in sperm 
(for c = 0 and x = 0.90 our model predicts a +/t heterozygotes 
paternity share of F + , = 0.28). Better estimates for this key pa- 
rameter of the model are clearly needed to further support the 
polyandry hypothesis. 

Female fitness advantage. +/t females in our popula- 
tion have a longer life expectancy, which is in marked contrast 
to what has previously been reported. An early study found simi- 
lar fitness advantages in both sexes (Dunn et al. 1958), but most 
recent publications show a different picture. In a detailed study 
on the fitness consequences of the t, Carroll et al. (2004) find fit- 
ness disadvantages for t carriers both in males and females (both 
in survivorship and reproductive output). They argue further that 
artificially high levels of competition in their enclosure popula- 
tions might amplify these differences. Indeed, we find indicators 
of high competition such as wounded adult males and females, 
and dead pups inside and outside nest boxes (presumably killed 
by conspecifics) in our study population — apparently not only 
among males, but surprisingly also among females. The fact that 
in contrast to the study of Carroll et al. (2004), individuals always 
had the opportunity to leave the population, makes this even more 
surprising. Reproductive skew turned out to be strong both in 
males and females (A. K. Lindholm, unpubl. data). Furthermore, 
more pups are produced each year than are able to recruit into 
the population (see Supporting information). However, we have 
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no information about what happens to mice that disperse out of 
the population. Our longevity estimates are based on nondispers- 
ing mice out of the barn, thus our result of genotype-dependent 
survival makes the assumption that there is no differential emi- 
gration. Even if our survival data are confounded with migration, 
it nevertheless seemed to be appropriate to use the longevity es- 
timates for the model, because they provide reliable information 
about the probability of an individual to establish in the population 
of interest. 

The conflict between different levels of 
selection. The implications of the combination of fitness 
effects between individual and mean population fitness are 
intriguing. Weak overall overdominance results if female positive 
and male neutral effects are merged together (see Fig. 4). It 
follows from this that maximum population fitness is reached 
if a considerable number of t haplotypes are still present in the 
population. This finding is counterintuitive, because distorter 
systems are usually characterized as "genomic parasites." In 
this case however, a certain amount of parasites is beneficial to 
the collective. This would suggest that there is an intermediate 



optimal rate of polyandry, as Arnqvist and Nilsson (2000) suggest 
in a meta-analysis on insects. Nevertheless, negative effects 
are still severe if two parasites occur in the same individual 
(tjt homozygosity). Meiotic drive makes these cases more 
frequent. The relation between equilibrium gene frequency and 
population fitness (Fig. 4) illustrates how nonrandom segregation 
can lead to outcomes that are no longer maximizing individual 
or population fitness. Sexual selection in the form of polyandry 
and sperm competition is a possible mechanism to recover losses 
in individual and population fitness by reducing the number of 
tjt zygotes produced. Hence, multiple mating can be seen as an 
adaptive strategy for the organism to increase individual fitness. 
The question however, if an internal t frequency equilibrium 
optimizing mean population fitness is reached, depends on the 
scenario. The observed t frequency dynamics do not offer a 
conclusive answer to this question, because t frequency at the 
end of the 5.5-year period falls between the population fitness 
optimum and t extinction. 

The conflict between levels of selection can also be seen 
as a conflict among the genes. Nonlinked genes can only spread 
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through populations when they increase individual fitness (as is 
the case if mean population fitness is maximized). In this respect, 
polyandry would represent an indirect solution to intragenomic 
conflict: by building up organisms that pursue a polyandrous mat- 
ing strategy, unlinked genes are able to prevent t haplotypes from 
reaching frequencies that are too high. However, the overall over- 
dominance effect is only weak and appears to be rather specific 
for the observed population. Thus, maximization of population 
fitness may not provide a general answer to the question of why t 
haplotypes generally occur at small but stable frequencies. 

The model predictions. Our specific case study provides 
a further example that the null model provided by Bruck ( 1 957) — 
even when allowing for drift — is not sufficient to explain the 
naturally observed frequencies. In contrast to the null model, our 
model shows that polyandry in the context of a meiotic drive 
system is a potential mechanism to explain the low t frequency 
paradox. 

There are plenty of reasons to assume that sperm competi- 
tion levels are substantially higher than assumed in scenario I. 
First, the assumption of a maximal number of two matings per 
female is clearly a simplification. Rolland et al. (2003), for ex- 
ample, showed that females can regularly change mating partners 
during one estrous cycle (13 out of 21 females received three 
ejaculations from two potential mating partners during one es- 
trous cycle). The probability p of having at least one / carrying 
male among all males mated increases with every mating partner. 
If k is the number of partners per estrous cycle and P+t(2) the 
proportion of heterozygote males in a population, this probabil- 
ity increases with k given p = (1 — P + , ( 2))*. This effect would 
increase the number of matings where the +/+ sperm compe- 
tition advantage plays a role. Second, we only looked at com- 
petitiveness effects dependent on the donor's genotype. This is of 
course a simplification, because other factors such as mating order 
usually have a major impact on reproductive success. The integra- 
tion of such factors would not change our predictions of t frequen- 
cies (assuming that the mating order remains random). However, 
it would strongly increase the discrepancy between multiple mat- 
ings and the actually measured multiple paternities. This stresses 
once more the urgency to estimate the specific factors determining 
male reproductive success in sperm competition in the future. 

Obviously, there is still the possibility that the forces included 
in the present model are insufficient to explain the observed t dy- 
namics. However, most of the potential factors that have been 
proposed in the last 50 years have been considered here. Popu- 
lation substructure appears rather unlikely given that migration 
within the barn is clearly too high to result in different demes. 
Even if there was a certain amount of substructure in the barn, 
Levin et al. (1969) showed that only small migration rates are 
needed to counterbalance the loss of t haplotypes by genetic drift 
as described in Lewontin and Dunn (1960). Inbreeding (Petras 



and Topping 1983) was excluded here as a possible factor (see 
Fig. S4). Migration of t carriers into the population is unlikely 
to affect t dynamics: We only found 5% unmarked adults (see 
Material and Methods) and genetic analyses did not identify im- 
migrants among them (A. K. Lindholm, unpubl. data). There is 
also the possibility that precopulatory female mating bias interacts 
with the here described postcopulatory effects of sperm compe- 
tition. At first sight, one might think that these two hypotheses 
could be exclusive. Why should females mate with multiple part- 
ners if they can recognize and discriminate between the different 
genotypes (as indeed suggested by Coopersmith and Lenington 
(1992))? However, female mating decisions might not always be 
as free and unconstrained as assumed. If the dominant male of 
their territory is a t carrier, females may not have any choice but 
to mate with him or risk infanticide (Perrigo et al. 1991). Results 
by Coopersmith and Lenington (1992) indicate that male domi- 
nance status is rated higher by females than genetic background. 
If one would assume that this preference is adaptive, this result 
would indicate that the costs for a female to mate exclusively 
with a subordinate male are even higher than the costs of mating 
with a +/t male. However, the correlation between male domi- 
nance and the t genotype is unclear. In arena experiments, +/t 
males tended to dominate +/+ males (Lenington et al. 1996), 
whereas Carroll et al. (2004) found +/+ males to be dominant in 
a seminatural context. In any case, an interaction of both multiple 
mating and precopulatory mate choice could definitely have more 
profound effects on equilibrium t frequency. It would be very 
interesting to investigate this further. For example, one might 
expect an increase in polyandry in the case of many t carrying 
males present in a population, as Price et al. (2008b) showed in 
D. pseudoobscura. 

General implications. Many have wondered why one 
cannot find modifiers in the genome that fight t distortion on 
a genetic level in naturally occurring house mouse populations 
(Ardlie 1998). The previously described X-linked SR system in 
D. pseudoobscura is one of the other cases where no modifiers 
have been found so far (Beckenbach 1996). The results presented 
here and those of Beckenbach (1996) and Ardlie (1998) sug- 
gest an alternative scenario: the possibility that there was never a 
need for such suppressors to evolve, because the t haplotype and 
the SR were already facing severe problems to be successful in 
polygynandrous mating systems. Thus, the mechanisms that such 
genetic elements use to make them strong intraejaculate competi- 
tors (reduced number of + sperm) can at the same time make them 
vulnerable to sperm competition, and hence, multiple mating. In- 
dividual behavior in the form of female mating decisions would 
in this case improve individual fitness (to the disadvantage of 
the t). Alternatively, the genes determining multiple mating could 
be regarded as another class of genetic modifiers of segregation 
distortion. 



EVOLUTION SEPTEMBER 2011 



2445 



ANDRI MANSER ET AL. 



The results presented here demonstrate that polyandry is 
a biologically plausible explanation for the low t frequency 
paradox. 
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